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Abstract 

We define a coupling of two baker maps tlirougli a 7r/2 rotation both in position and in momen- 
tum. Tlie classical trajectories thus exhibit spiraling, or loxodromic motion, which is only possible 
for conservative maps of at least two degrees of freedom. This loxodromic baker map is still hyper- 
bolic, that is, fully chaotic. Quantization of this map follows on similar lines to other generalized 
baker maps. It is found that the eigenvalue spectrum for quantum loxodromic baker map is far 
removed from those of the canonical random matrix ensembles. An investigation of the symmetries 
of the loxodromic baker map reveals the cause of this deviation from the Bohigas-Giannoni-Schmit 
conjecture. 

PACS numbers: 
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I. INTRODUCTION 



The essential feature of classically chaotic motion is hyperbolicity, that is both stretching 
and squeezing, within a confined region of phase space. The simplest textbook examples, as 
far as (conservative) Hamiltonian systems are concerned, are two-dimensional {2D) maps, 
which may be pictured as Poincare sections of motion in a AD phase space. However, hy- 
perbolicity is in no way confined to lower dimensions. In particular, a conservative return 
map in four dimensions may exhibit spiralling, or loxodromic motion onto a fixed point. 
This would not be allowed in two dimensions since it would violate area conservation, but is 
compensated in four dimensions by outward spiralling in a complementary plane. The lin- 
earization of hyperbolic conservative maps of arbitrarity high dimension can be decomposed 
into products of simple hyperbolic motion in eigenplanes together with loxodromic motion 
in 4i5-subspaces 

Hyperbolic motion is generally desribed by horseshoe maps in the qualitative theory of 
dynamical systems . The idea is to describe by symbolic dynamics the successive visits 
of a trajectory to a well chosen finite set of domains of phase space. The problem, as far as 
conservative systems are concerned, is that most orbits are not restricted to this finite set 
of domains, so that the horseshoe map itself is not conservative. This is the motivation for 
the develoment of a specially convenient model known as the baker map. 

The baker map displays the essential features of classically chaotic motion in such a 
simplified form as to be almost a caricature 0, S]. It may be described as a space-filling 
horseshoe map that linearizes the hyperbolic motion. Thus, essentialy chaotic motion is 
exhibited, which can be followed through very simple computations. The binary symbolic 
dynamics propagates vertical strips in the unit square onto horizontal strips. The vertical 
and horizontal rectangles become narrower, for longer binary codes pertaining to multiple 
iterations of the baker map. The primary digit specifies either of the two half-squares, so 
that it is responsible for a coarse-grained description of the motion. The way in which 
horizontal strips are piled up by the mapping is also determined by the primary binary digit 
(see Fig. OP). 

Higher dimensional generalizations of the baker map have recently been presented joj. 
The products of horizontal rectangles for each degree of freedom define hyper-parallelepipeds 
that are stretched and squeezed into horizontal parallelepipeds. Even though the increased 
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dimension admits a much richer range of possibihties for restacking these domains, none of 
these leads to spiraUing motion. Here we show how a rotation can be added which leads to 
spirals, without violating any of the essential attributes of generalized baker maps. 

The various quantization schemes that have been proposed for the baker map generally 
respect the main division of the square into a pair of vertical rectangles, which split the posi- 
tion states into Hilbert subspaces. These are mapped respectively onto a pair of momentum 
state subspaces according to the primary binary digit. The discrete and finite nature of the 
Hilbert space prevents the association of longer codes to ever thinner strips Thus, to 
some extent, the primary digit is even more important in quantum mechanics. 

Initially, it may have seemed that the introduction of binary symbols to describe evolving 
quantum systems could only serve as an artificial prop for the study of the semiclassical limit. 
However, the growing interest in quantum computation has brought the binary structure to 
the fore. In this context, the quantum baker map has reemerged as an ideal simplified model. 
Schack and Caves Q| have indeed developed an entire class of quantum baker maps, based 
on the 2^-dimensional Hilbert space of_N qubits, which include the original quantizations 
of Balazs and Voros [9] and Saraceno jlO] as special cases. The Schack-Caves family was 
recently further enlarged by Ermann and Saraceno [11], who also explained the general 
structure of all these maps. We here adopt the Balazs- Voros-Saraceno 10| map throughout. 

It is important at this stage to distinguish two alternative forms of classical-quantum 
correspondence for individual baker maps. In the original quantizations ol, 1^, the classical 
phase space is two-dimensional, a torus. It is the chaotic motion in this finite region that 
gives rise to the repulsion of the quantum quasi-energy levels, characteristic of the spectra 
of random matrices, according to the Bohigas-Giannoni-Schmit (BGS) Conjecture ||1^]. If 
the quantization is chosen with a value of Planck's constant such that the Hilbert space 
has exactly 2^ states, then it is possible to reinterpret the quantum system to be the 
tensor product of two-level systems. This would correspond to a hypercube in a 2N- 
dimensional phase space, but it would be stretching it to ascribe a classical correspondence 
to this quintessentially quantum system. We shall here keep to the original quantization 
schemes and, hence, interpret each quantum baker map as corresponding to a classical 
two-dimensional phase space, so that the coupling of two baker maps corresponds to a 
four-dimensional phase space. The four possible combinations of primary binary digits, 
(0, 0), (0, 1), (1, 0), (1, 1), then correspond to four squares in the classical position space. 
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(gi, ^2) and hence four parallelepipeds in the unit 4-cube, (gi, q2,Pi,P2) in which the classical 
motion is defined. 

The interaction between baker maps here involves a rotation of the four parallelepipeds 
into which the primary digits of both maps divide the four dimensional phase space. Choos- 
ing this to be of 7r/2, merely exchanges the parallelepipeds, so that the doubled binary 
decription is respected. The interesting point is that the equilibrium and all the classical 
periodic orbits become "loxodromic" , i.e. the positions spiral outwards, while the momenta 
spiral inwards, because the eigenvalues of the stability matrix are general complex numbers: 
Neither are they real, nor do they lie on the unit circle. Obviously, this example models more 
closely a two dimensional motion for a single particle than any easily conceivable dynamics 
for a pair of one-dimensional particles. Loxodromic behaviour does not arise for equilibria 
of simple p"^ + V Hamiltonians, but it is a generic possibility for the Poincare map around 
periodic orbits in higher dimension. So far, not much effort has been made towards quan- 
tizing loxodromic motion, except for higher dimensional cat maps ISj, another simplified 
model. 

The analysis of the loxodromic baker map is preceeded by a brief review of the ordinary 
two-dimensional baker map and its quantization. We point out the possibilities for different 
pilings that lead to the simplest couplings of baker maps [6|]. Section 3 then presents the 
loxodromic baker map and its quantization. It is shown in section 4 that, contrary to the 
generalized baker maps in j^, the quantum eigenvalue spectrum does not obey a generic 
random matrix ensemble. The reason for this failure of the BGS conjecture is shown to 
be that the 7t/2 rotation commutes with the component baker map. This leads to a non- 
trivial decomposition into a tensor product of two maps, so that the quantum quasi-energy 
spectrum in this case is highly degenerate. 



II. REVIEW OF CLASSICAL AND QUANTUM BAKER MAPS 

In this section we present the well known classical and quantum ingredients of baker 
maps. The classical baker transformation is an area preserving, piecewise-linear map, b : 
(Po^Qo) ^ (Pi^Qi) of the unit square (periodic boundary conditions are assumed) defined as 

Pi = ^(Po + eo), gi = 2go-eo; (1) 
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FIG. 1: The baker map propagates a pair of vertical rectangles onto a pair of horizontal rectangles. 
Inside each separate region the evolution is linear. 

where eo = [2go], the integer part of 2go- The two regions into which the phase space is 
separated and their evolution is shown in Fig. [H 

This map is known to be uniformly hyperbolic, the stability exponent for orbits of period 
L being Llog2. Moreover it admits a useful description in terms of a complete symbolic 
dynamics. A one to one correspondence between phase space coordinates and binary se- 
quences, 

(p, g) ^ . . .e_2e-i ■ eoeie2 . . . , ei = 0, 1 , (2) 

can be constructed in such a way that the action of the map is conjugated to a shift map, 
i.e., to shifting the dot to the right in the binary sequence above. The symbols are assigned 
as follows: ej is set to zero (one) when the i-th iteration of (p, q) falls to the left (right) of 
the line g = 1/2, i.e. [2gj] = e^. Reciprocally, given an itinerary 

. . .e_2e_i ■ eociea . . . , (3) 

the related phase point is obtained through the specially simple binary expansions 

oo ^ oo g 

9 = E^' p = E^- (4) 

i=0 ^ i=l ^ 

Once the dynamics has been mapped to a shift on binary sequences it is very easy to analize 
the dynamical features of the map. In particular, periodic points are associated to infinite 
repetitions of finite sequences of symbols. 

Due to its piecewise linear nature, the baker map admits a (mixed) generating function 
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which is a piecewise bihnear form, 

Weoipi, go) = 2pigo - eoPi - eogo , eo = 0, 1 . (5) 
It is not defined on the whole space pi-qo but on the classically allowed domains 

i?o = [0, 1/2] ® [0, 1/2] and i?i = [1/2, 1] ® [1/2, 1] . (6) 

Though the above generating function will be the starting point for quantization, it must 
be remembered that it only provides an implicit formula for each iteration of the classical 
baker map. 

With respect to the quantum map, we will follow the original quantization of Balazs 
and Voros 9], as later modified by Saraceno |lOj] to preserve in the quantum map all the 
symmetries of its classical counterpart. In the mixed representation the baker's propagator 
can be written as a D x D block-matrix [D even): 

/ 



{Pni\BD\qn) 



Gd/2 



V 



G 



(7) 



D/2 



where position and momentum eigenvalues run on a discrete mesh with step 1/D = h {h 
Planck's constant), so that 



qn = {n + l/2)/D , pm = {m + l/2)/D , 0<n,m<D-l ■ 
and is the antiperiodic Fourier matrix, which transforms from the q to the p basis. 



(8) 



Gd = {Pm\qn) = (l/V:D)e-2-^^'"^" 



(9) 



It will be useful to consider this matrix as the g- representation of the Fourier operator Gd. 
The propagator for the baker map has the standard structure of quantized linear sym- 



plectic maps 



{pm\BD\qn) 



otherwise . 



(10) 



In this quantization, only those transitions are allowed that respect the rule [2pm] = [2g„], a 
refiection of the classical shift property. To be able to iterate the quantum baker map, the 



state must be brought back to the position representation. This is achieved by an inverse 
Fourier transform, so that the rnatrix ([7]) is multiphed by G^. 

Following Schack and Caves |8'], we can reinterpret this quantum map as the evolution of 
N qubits if the dimension of the Hilbert space satisfies D = 2^ . Then the position states 
can be defined as product states for the qubits in the basis, 

\qn) = |ei) ® |e2) ® ...|eAr), 
where n has the binary expansion 

N 

n = ei...e„ = ^e,2^-^' (12) 
i=i 

and Qn = {n + 1/2) / D = 0-ei...eArl. The connection with the classical baker map is specified 
by the symbolic djTiamics. The bi-infinite strings ([3]) that determine points in the unit square 
are made to correspond to sets of orthogonal quantum states. Half of the position states lie 
in either of the two rectangles, Rq or Ri defined in (jH]), which correspond respectively to 0, 
or 1 eigenstates of the principal qubit. 

The full unitary operator for the quantum baker can be written explicitly as 

BD = Go'[h®Gn/2], (13) 

where I2 is the unit operator for the first bit and G0/2 is the Fourier operator on the re- 
maining qubits. The operator in the square brackets preserves the first qubit, while evolving 
separately the remaining qubits, within each domain R^^. It is the final inverse Fourier op- 
erator that mixes the principal qubit in with the rest, because it acts globaly on the states 
in both domains. This step is not explicitated in the mixed representation ([7j). 

So far we have only allowed for a single possibility in which to stack the rectangles in 
the baker transformation, but alternative to ([1]), the classical map, b' : {po,qo) — > {PiyQi), 
shown in Fig. [2], 

Pi = ^(Po + 1 - eo) , gi = 2go-eo, (14) 

has very similar properties. This variation is tantamount to reversing the primary classical 
bit and it corresponds to the quantum map 13', represented by the matrix: 

{Pm\B'jj\qn) =1 . (15) 

\Gd/2 ) 
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FIG. 2: The alternative baker map also propagates a pair of vertical rectangles onto a pair of 
horizontal rectangles, but their stacking is interchanged. 

In other words, we here substitute the operator I2, which acted on the principal qubit, by 
X2 represented by the Pauli matrix: 



so that % = G^^[X2 0Gd/2]- 

Of course, we are free to substitute any other unitary operator acting on the primary 
qubit, but it is only B' that can be interpreted classically as an equivalent alternative piling 
of the baker map. We can still split up the evolution into domains equivalent to (Q, though 
with a different matching of go and pi segments. The classical evolution within each domain 
is again determined by classical generating functions like which becomes the exponent 
of the propagator. 

If we allow a finite probability for a stacking fault in the classical baker map, the evolution 

acquires a random component. However, if the piling order depends on the coarse-grained 

position of another baker map, the overall motion will again be purely deterministic. Being 

that each baker map is thoroughly chaotic, it will be hard to distinguish the random motion 

I ^ 

of one of the components taken on its own from a truly stochastic system. In [6] we allowed 
the baker to choose between these alternatives, depending on its interaction with another 
baker. Though this led to surprisingly rich structures, no spiraling motion could be generated 
in the corresponding classical map. This possibility is now realized by rotating the domains 
into which the baker's phase space is partitioned. 




(16) 
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FIG. 3: The partition of phase space into disjoint domains is respected by the product of baker 
maps. 4D-parallelepipeds, projecting onto squares in the position space are transformed onto 
horizontal parahelepipeds that project onto squares in momentum space. 

III. DEFINITION OF THE LOXODROMIC BAKER MAP 



The direct product of a pair of classical baker maps, ® b^, is itself a generalized baker 
map defined within a 4D-hypercube, {p^, q^) (p^, q"^)- Indeed the motion is hyperbolic and 
the partition of phase space into disjoint domains is again respected. These are specified by 
a doubled binary code, eg = cg^ ® eo^? defining ^;ertzca/ 4D-parallelepipeds that project onto 
squares in the position space (g^, g^), while covering the full available momentum space. The 
effect of a single iteration of this product map is to transform these domains onto horizontal 
parallelepipeds that project onto squares in momentum space, as shown in Fig. [3l 

The main difference between this product map and the generalized bakers defined in jc] 
is that in the latter case the symbolic mapping is no longer the identity. One could certainly 
define a generalized baker by repiling the product maps within each domain so as to obtain 
the symbolic evolution displayed in Fig. HI though this possibility was not singled out in . 

However, the same effect (as far as the symbolic evolution is concerned) can be obtained 
by combining the product baker with a canonical point transformation, r™2 5 that simply 
rotates both the position and the momentum spaces around the midpoint of the hypercube, 
m = (1/2, 1/2, 1/2, 1/2), as shown in Fig. B 

This loxodromic baker map, ° [b^ ® b^] = bL : (po^Qo) ~^ (Pi^Qi) is given explicitely 
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FIG. 4: Possible repiling of the partition of the domains of the product of baker maps. 





FIG. 5: A canonical point transformation, r^g; that simply rotates both the position and the 
momentum spaces around the midpoint of the hypercube, m = (1/2, 1/2, 1/2, 1/2), has the same 
effect on the partition of the product of baker maps as the repiling in Fig. [H 
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FIG. 6: Typical classical trajectories of the loxodromic baker map spiral out of the fixed points in 
the q-plane and into the fixed points in the p-plane. Every four iterations they return to the same 
radial line, because the rotation angle is 7r/2. 

This is no longer a product map within each separate domain and Fig. [6] shows how an 
orbit spirals onto a fixed point in momentum space and spirals outwards in position space. 

The loxodromic baker map can also be described by symbolic dynamics. This now has a 
four letter alphabet, two binary digits, instead of one for the ordinary baker map. Further- 
more, the relation between the symbolic description and the phase space trajectories is not 
so direct. The details may be obtained in jl6|]. 

The quantum version of the loxodromic baker map is defined in strict correspondence with 
the foregoing classical construction. This is always possible for linear classical evolutions, i.e. 
we can define a quantum unitary transformation, -R^2 5 which corresponds to the classical 
transformation, r^M, so that, in combination with the definition of the usual baker map 
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(fT3|) . we then stipulate that 

Bl = RT/2 BD®Bn. (21) 

The only thing that is now missing is the full definition of the rotation operator in the 
above formula. Recalling that the classical rotation is about the midpoint m in phase space, 
it follows that, if the rotation operator about the origin by an angle a is 

= exp [-ia{p2qi - q2Pi)] , (22) 

then 

= Tm R-a T-m, (23) 

where is the quantum operator corresponding to a uniform translation (or displacement) 
by the phase space vector ^ (see e.g.fisl). 

A considerable simplification is possible for the 7i/2 rotation, which is the only one allowed 
for a generalized baker map. In this case, it is easy to verify that the classical rotation about 
the midpoint, m, can be written as a composition 

= t, o r°/2, (24) 

where is the classical translation by the vector r] = (1,0, 1,0) and r°y2 the tc/2 rotation 
around the origin. Therefore, we can quantize this form directly so as to obtain 

= Tr, Rl/2- (25) 



IV. SYMMETRIES AND LEVEL REPULSION 



The celebrated BGS conjecture predicts the correspondence between chaotic classical mo- 
tion and quantum level repulsion. Specifically, for quantum maps, the eigenvalue spectrum 
of the unitary operator should have fiuctuations in agreement with the canonical CUE, or 
COE ensembles of random matrix theory (RMT). Exceptions do arise because of symme- 
tries, or even more subtle arithmetic anomalies which may be hidden in the definition of 
the map. However, the BGS conjecture has become so well established that its negation is 
taken as a sure sign that hidden structure should be sought out. 

In relation to the baker map, it should be noted that the original quantization by Balazs 
and Voros did not conform properly to the conjecture. Saraceno then discovered that 
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FIG. 7: Distribution of level spacings -normalized to unit average- for the repiled baker maps 
(full dimension equal to 512; big dots). We also display the predictions of random matrix theory 
(dotted line: COE, dashed: CUE) and the Poisson distribution (full line). 

this form of quantization did not incorporate an exact classical symmetry and proposed a 



modification lOj. The individual symmetry classes of this new map are satisfactory examples 



of the BGS conjecture as can be verified in 15|. 



Even though the direct product of two baker maps is certainly a chaotic classical sys- 
tem, its quantization does not satisfy the BGS conjecture. Indeed, the eigenvalues will 
be the products of the baker eigenvalues, leading to degeneracies instead of level repul- 
sion. In general, the addition of some interaction between the two maps, while taking the 
precaution to preserve the hyperbolicity of the classical motion, should lead again to an 
RMT spectrum. In the case of the generalized "repiled" baker maps studied in 6[ the level 
spacing distribution, shown in Fig. [71 does not comply with the BGS conjecture. Its in- 
termediate nature between a Poisson spectrum and a COE spectrum suggests the nearness 
of a symmetry, as was found for the original Balasz-Voros quantization of the baker map. 
However, Fig. [8] clearly rules this possibility out in the case of the loxodromic baker map. 
The Poisson profile indicates independent samples and, furthermore, there are many degen- 
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FIG. 8: Histogram of level spacings for eigenangles of a loxodromic baker map with dimension 
equal to 1296. The smooth curve corresponds to the Poisson spectrum. 

eracies, just as with the simple product map. Therefore we must seek out the symmetry 
that is responsible. It turns out that this involves the reflection symmetry of the classical 
3aker map, R : {q,p) R{q,p) = (1 — g, 1 — p), included in the Saraceno quantization 
10(1 . Of course, both factors of the product map, ® b^, commute with the respective 
reflections, Ri and R2. Another obvious symmetry of the product map is the permutation, 
P : g2,P2) — ^ P{(li-,Pii<l2-,P2) = iq2,P2,qi,Pi)- Combining these, we then obtain: 

RiP(gi,Pi,g2,P2) = Riiq2,P2,quPi) = (1 - q2, i -P2,qi,pi) = r^/2(gi,Pi,g2,P2)- (26) 

Thus, we flnd that the special nature of the 7r/2 rotation, that was chosen because it respects 
the geometrical divisions of the baker map, can now be reinterpreted as a particular product 
of the symmetries of b^ (8) b^. 

It is then neither surprising, nor hard to prove, that itself is a symmetry of the 
classical and the quantum product baker map. The consequence is that i?^/2 and i?D ® Bd 
have common eigenfunctions and they must share these with B^, as defined by (|2T1) . Thus, 
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if we divide the spectrum of the product map into the symmetry classes of the quantum 
rotation R-„/2, with eigenvalues {+1, —1, +i, — the eigenvalues of Bl in each class are just 
the corresponding multiple of the eigenvalues of the quantized product map. In conclusion, 
the loxodromic baker map belongs to the class of quantized chaotic maps which fail to satisfy 
the BGS conjecture because of a symmetry. 

V. DISCUSSION 

The quantization of spiraling chaotic motion has not received the attention it deserves, 
because, so far, work has concentrated on two-dimensional maps. For systems with many 
degrees of freedom, it must be borne in mind that Poincare maps will generally exhibit 
spiraling behaviour along with the more familiar features studied in two-dimensional maps. 
The complexity of higher dimensional motion encourages the search for specially simple 
models with which to study the classical-quantum correspondence and the semiclassical limit. 
The generalized baker map presented here was an obvious candidate, but the restriction to 
rotations of 7r/2, amounting roughly to a permutation of the variables between the two 
coupled maps, turns out to be too severe. The commutation of this rotation with the baker 
product map leads to a totally unsatisfactory eigenvalue spectrum for a chaotic quantum 
map. This is in sharp contrast to even simpler couplings developed for baker maps in 

Another possibility that has also been investigated are loxodromic cat maps, i.e. linear 
maps on the 4-torus There is then a richer range of possibilities, but still there are 

restrictions to the distribution of eigenvalues common to all cat maps. These arise from 
the necessity of quantum cat maps to be periodic, albeit the period may be very high. 
Presumably, the search for simple loxodromic quantum models will only heat up when a 
loxodromic algorithm is invented for quantum information theory. 
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